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ABSTRACT 


KROLL,  JAMES  T.  A  Case  Study  of  a  Mid-Atlantic  Coastal 
Front .( Under  the  direction  of  DAVID  BARBER  and  GERALD 
WATSON ) . 

^A  unique  case  of  mid-Atlantic  coastal  f rontoqene si s  on 
December  21-22,  1983  is  examined  to  determine  the  evolution 
of  the  three  dimensional  structure  of  this  mesoscale 
phenomena.  Analysis  of  synoptic  and  mesoscale  surface  maps 
suqqest  a  coastal  front  which  moves  steadily  inland. 

However,  review  of  temperature  and  wind  profiles  at 
individual  stations  reveals  that  f rontoqenesi s  actually 
occurs  over  the  land  in  North  Carolina  and  southern 
Virqinia.  Laqranqian  and  Eulerian  f rontoqenesi s 
calculations  confirm  this  findinq.  The  f rontoqenesi s 
process  is  found  to  be  larqely  a  result  of  the  aqeostrophic 
flow  within  the  cold  wedqe  which  results  in  weak  cold  air 
advection  over  the  land.  This  northerly  flow  combines  with 
a  stronq  southerly  flow  over  the  Atlantic  to  increase  the 
temperature  qradient  near  the  coastline. 

Analysis  of  kinematic  variables  indicates  reqions  of 
Lelatively  stronq  converqence  and  cyclonic  vorticity  in  the 
reqion  of  the  coastal  front.  Geostrophic  deformation  maxima 


coincide  with  the  l  idqe  of  the  anticyclonic  bulge  east  of 
the  Appalachians.  Three-dimensional  Laqranqian 
f rontoqenesis  computations  identify  the  "shear"  and 
"confluence"  terms  as  f rontoqenetical  within  the  frontal 
zone.  The  "tilting"  term  is  f rontol yt ica 1  within  the 
frontal  zone.  Precipitation  analyses  identify  the  cold  air 
side  of  the  coastal  front  as  favorable  for  enhanced 
precipitation.  Relationships  between  preciptation  bands  and 
convective  instability  are  also  discussed.^ 

To  test  the  hypothesis  that  the  inland  ageostrophic 
flow  is  a  result  of  problems  with  the  reduction  of  pressures 
to  sea  level,  a  terrain  level  qeostrophic  wind  is  used  to 
predict  wind  flow.  The  terrain  level  winds  failed  to 
produce  a  better  estimate  of  actual  winds  than  that  of  the 
sea  level  qeostrophic  winds. 
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ABSTRACT 


KROLL,  JAMES  T.  A  Case  Study  of  a  Mid-Atlantic  Coastal 
Front. (Under  the  direction  of  DAVID  BARBER  and  GERALD 
WATSON ) . 

A  unicue  case  of  mid-Atlantic  coastal  f rontogenesi s  on 
December  21-22,  1983  is  examined  to  determine  the  evolution 
of  the  three  dimensional  structure  of  this  mesoscale 
phenomena.  Analysis  of  synoptic  and  mesoscale  surface  maps 
suqqest  a  coastal  front  which  moves  steadily  inland. 

However,  review  of  temperature  and  wind  profiles  at 
individual  stations  reveals  that  f rontogenesis  actually 
occurs  over  the  land  in  North  Carolina  and  southern 
Virqinia.  Laqranqian  and  Eulerian  f rontoqenesis 
calculations  confirm  this  findinq.  The  f rontoqenesis 
process  is  found  to  be  larqely  a  result  of  the  ageostrophic 
flow  within  the  cold  wedge  which  results  in  weak  cold  air 
advection  over  the  land.  This  northerly  flow  combines  with 
a  strong  southerly  flow  over  the  Atlantic  to  increase  the 
temperature  gradient  near  the  coastline. 

Analysis  of  kinematic  variables  indicates  regions  of 
relatively  strong  convergence  and  cyclonic  vorticity  in  the 
region  of  the  coastal  front.  Geostrophic  deformation  maxima 
coincide  with  the  ridge  of  the  anticyclonic  bulge  east  of 
the  Appalachians.  Three-dimensional  Lagranqian 
f rontogenesis  computations  identify  the  "shear"  and 


"confluence”  terms  as  f rontogenetical  within  the  frontal 
zone.  The  "tilting"  term  is  f rontolyt ical  within  the 
frontal  zone.  Precipitation  analyses  identify  the  cold  air 
side  of  the  coastal  front  as  favorable  for  enhanced 
precipitation.  Relationships  between  preciptation  bands  and 
convective  instability  are  also  discussed. 

To  test  the  hypothesis  that  the  inland  ageostrophic 
flow  is  a  result  of  problems  with  the  reduction  of  pressures 
to  sea  level,  a  terrain  level  geostrophic  wind  is  used  to 
predict  wind  flow.  The  terrain  level  winds  failed  to 
produce  a  better  estimate  of  actual  winds  than  that  of  the 
sea  level  geostrophic  winds. 
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INTRODUCTION 


1.1  General  Comments 

In  the  colder  months  of  the  year,  the  synoptic  pattern 
over  the  Eastern  United  States  frequently  enhances  the 
development  of  two  mesoscale  phenomena,  namely,  the  cold 
wedge  and  the  coastal  front.  The  cold  wedge,  also  known  as 
the  "Appalachian  damming  effect,"  is  the  damming  or  trapping 
of  a  shallow  layer  of  relatively  cold  air  east  of  the 
Appalachian  mountains.  The  coastal  front  is  an  intense 
mesoscale  baroclinic  zone  that  develops  over  the  relatively 
warm  adjacent  coastal  waters. 

The  synoptic  pattern  responsible  for  these  mesoscale 
events  is  dominated  by  a  cold  anticyclone  centered  over  New 
England  or  southeastern  Canada.  A  striking  feature  which 
develops  within  the  cold  wedge  is  the  maintenance  of 
northerly  flow  over  the  inland  region,  east  of  the  mountain 
range.  This  flow  is  contrary  to  what  is  expected  by  the 
analyzed  pressure  pattern.  This  ageostrophy  is  especially 
pronounced  along  the  base  of  the  mountain  range  where 
surface  wind  direction  occasionally  differ  from  the 
geostrophic  direction  by  more  than  90  degrees. 

Previous  studies  indicate  that  the  ageostrophic  flow 
within  the  cold  wedge  appears  to  be  a  significant  factor  in 
developing  and  maintaining  the  coastal  front.  While 
northeasterly  flow  dominates  inland,  the  flow  over  the  ocean 
gradually  shifts  to  the  southeast.  The  coastal  front 
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develops  in  the  boundary  zone  between  the  opposing  flows. 
Frontal  development  is  enhanced  as  the  air  over  the  Atlantic 
traverses  long  fetches  of  warmer  ocean  waters  causing 
temperature  modification  by  sensible  heating.  This 
ultimatly  results  is  an  increase  of  the  overall 
onshore-offshore  temperature  gradient. 

The  occurrence  of  these  phenomena  is  significant  to  the 
synoptician.  Freezing  rain  is  frequently  associated  with 
these  events  during  the  winter  months.  The  coastal  front 
itself  is  often  the  dividing  line  between  frozen  and  liquid 
precipitation.  Considering  that  the  frontal  temperature 
gradient  may  be  as  great  as  1  K/Km  over  a  20-30  Km  wide 
zone,  the  prediction  of  frozen  versus  non-frozen 
precipitation  regions  is  difficult.  Adding  to  the  dilemma 
is  the  inability  of  the  operational  forecast  models  to 
handle  correctly  this  boundary  layer  phenomena.  The 
forecast  models  tend  to  eliminate  the  high  pressure  wedge 
and  establish  a  warming  trend  over  the  coastal  reqions. 

This  rarely  occurs. 

Bosart(1972)  identified  New  England  and  the  Carolinas 
as  two  geographical  regions  with  a  marked  preference  for 
coastal  front  development.  Bosart  and  Korty(1976)  linked 
the  Carolina  coastal  front  with  the  redevelopment  of 
cyclones  moving  norhteastward  from  the  Gulf  of  Mexico.  To 
date,  however,  detailed  analyses  of  coastal  fronts  are 
limited  to  one  analysis  in  each  of  these  preferred  regions. 


1.2  Objectives 

As  a  preliminary  study  of  the  East  Coast  Cyclone 
Project  (GALE),  this  research  presents  a  case  study  to 
illustrate  a  mid-Atlantic  coastal  front  event.  Efforts 
concentrate  on  depicting  the  surface  and  three-dimensional 
structure  of  the  coastal  front  as  well  as  identifying  the 
processes  which  affect  frontal  evolution.  Specifically,  an 
objective  analysis  scheme  is  used  to  derive  pertinent 
meteorological  parameters  (i.e.  divergence,  gradients, 
advections,  etc.).  Eulerian  and  Lagrangian  forms  of  the 
f rontogenesis  equation  are  also  applied.  Relationships 
between  these  results  and  frontal  evolution  and  movement  are 
established . 

The  relationships  between  frontal  structure  and  weather 
distibution  are  also  examined.  Potential  instability 
analyses  are  examined  in  conjunction  with  the  precipitation 
distribution.  In  addition,  terrain-level  geostrophic  winds 
are  compared  with  geostrophic  winds  based  upon  pressures 
reduced  to  sea  level  as  an  alternative  to  predicting  wind 
flow  within  the  cold  wedge. 

The  selected  case  is  not  associated  with  any  marked 
cyclogenesis.  The  event  is  a  moderately  intense  coastal 
front  which  passes  through  reporting  stations  80-120  Km 
inland.  Because  the  frontal  zone  passes  through  a  dense 
network  of  surface  reporting  stations,  a  more  complete 
analysis  of  frontal  evolution  may  be  obtained. 


1.3  Literature  Review 
1.3.1  General  Frontal  Studies 

J.  Bjerknes( 1919 )  is  generally  credited  with 
documenting  the  first  significant  description  of  surface 
front  structures  and  the  relationship  to  extratropical 
cyclones.  He  noted  two  lines  of  convergence  intersecting  at 
the  cyclone  center  which  bounded  the  warm  air  region. 

Bjerkes  and  Solberg ( 1922 )  elaborated  upon  the  cyclone  model 
to  include  its  entire  life  cycle.  According  to  the  theory, 
the  polar  front  was  a  zero-order  discontinuity  in  density  or 
temperature  which  intersected  the  earth's  surface. 

With  the  advent  of  upper  air  soundings,  Bjerknes( 1932 ) 
discovered  first-order  discontinuities  in  the  upper  level 
frontal  boundaries.  Instead  of  sharp  contrasts  in 
temperature,  he  discovered  a  sharp  contrast  in  temperature 
gradients.  This  led  to  the  concept  of  fronts  as  zones  of 
transition  and  not  lines  of  abrupt  change. 

Petterssen ( 1940  )  extensively  discussed  f rontogenesis  in 
relationship  to  wind  fields.  The  concept  of  f rontoqenesis 
involves  the  juxtaposition  of  two  differing  air  masses 
leading  to  the  creation  of  a  discontinuity  in  the 
temperature  field.  Defining  " frontogenesis"  as 

F  =  (d/dt) 1  VS  I  ,  1.1 

where  VS  is  the  magnitude  of  the  gradient  of  a 


conservative  property  S,  f rontogenesis  occurs  where  the 
gradient  of  S  increases  most  rapidly.  Petterssen  displays 
deformative,  divergent  and  rotational  fields  in  relation  to 
f rontogenesis.  He  concludes  that  deformative  fields  are 
f rontogenetical  when  the  isotherms  lie  within  a  45  degree 
angle  of  the  axis  of  dilatation.  Isotherms  lying  at  a 
larger  angle  to  this  axis  make  the  deformative  field 
f rontolytical .  In  addition,  he  resolved  the  inherent 
f rontogenetical  property  of  a  field  of  convergence  and 
determined  that  pure  rotational  fields  were  neither 
f rontogenetical  nor  f rontolytical . 

Prior  to  this  time,  f rontogenesi s  concepts  were  limited 
to  the  horizontal  plane.  Vertical  motions  in  the  vicinity 
of  surface  fronts,  except  for  forced  upslope  or  downslope, 
were  assumed  to  be  negligible.  The  expanding  network  of 
upper-level  soundings  provided  the  opportunity  forfrontal 
analysis  at  upper  levels  as  well  as  derived  values  of 
vertical  motion  in  the  atmosphere.  Miller (1948)  altered  the 
concepts  of  frontogenesis  when  he  redefined  it  as 

F  =  d  S VS I /dt  1.2 

In  addition  to  extending  the  f rontogenetical  concepts  into 
three  dimensions,  Miller's  equation  also  viewed 
frontogenesis  from  the  parcel  following  perspective.  This 
perspective  has  dominated  frontal  analysis  literature  since 
that  time  . 

Palmen(1948)  illustrated  the  temperature  and  wind 
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distributions  in  the  vicinity  of  the  polar  front.  The  polar 
front  was  identified  as  a  region  of  strong  stability  with  an 
associated  wind  maximum  at  upper-levels.  An  extension  of 
this  work  by  Palmen  and  Newton (1948)  verified  the  existence 
of  a  jet  stream  above  the  frontal  layer  in  the  vicinity  of 
the  tropopause.  The  polar  front  itself  had  a  varying  slope 
of  1/120  at  500mb  compared  to  1/175  near  the  earth's 
surface.  In  addition,  the  mean  horizontal  wind  shear  was 
cyclonic  throughout  the  entire  frontal  zone  with  a  distinct 
region  of  divergence  at  the  southern  edge  of  the  jet  stream. 

Miller's  concept  of  three-dimensional  frontogenesis  was 
applied  at  the  500mb  level  by  Reed  and  Sanders( 1953 )  and 
again  by  Reed (1955).  Results  indicated  that  the  horizontal 
gradient  of  vertical  motions  (referred  to  as  the  "tilting" 
term)  was  primarily  responsible  for  frontogenesis  at  that 
level.  Strong  upper-level  subsidence  at  the  warm  air  side 
of  the  front  was  largely  responsible  for  the  dominance  of 
that  term. 

Applying  the  same  concepts,  Newton (1954)  evaluated 
factors  affecting  frontogenesis  at  multiple  levels  in  the 
atmosphere.  He  demonstrated  that  factors  affecting 
frontogenesis  varied  at  different  atmospheric  levels.  The 
"tilting"  term  was  a  significant  f rontogenetic  factor  at 
mid-levels  as  was  determined  by  Reed  and  Sanders.  On  the 
contrary,  divergence  and  thermal  advection  within 
deformation  fields  dominated  f rontogenetic  effects  at  the 
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surface  and  at  the  level  of  maximum  winds.  Newton  suggested 
that  vertical  motions  near  the  surface  are  damped  enough 
such  that  their  effect  on  surface  f rontogenesis  was 
negligible.  Also/  surface  f rontogenesi s  was  usually 
associated  with  increased  vorticity  and  vorticity  changes 
were  affected  by  divergence  and  friction.  These  concepts 
had  been  hypothesized  earlier  by  Petterssen  and 
Austin ( 1942 ) . 

The  processes  of  intense  surface  f rontogenesis  were 
investigated  by  Sanders (  1955  ) .  Vertical  cross  sections  of 
the  "confluence"  and  "tilting"  terms  of  Miller's 
f rontogenesis  eguations  were  developed.  These  patterns 
indicated  that  the  "confluence"  term  was  highly 
f rontogenetical  within  the  frontal  zone  near  the  surface 
(probably  due  to  the  extreme  confluence  and  temperature 
gradients  occurring  at  the  surface).  On  the  other  hand,  the 
"tilting"  term  displayed  high  f rontogenetical  values  in  the 
warm  air  adjacent  to  the  front  while  the  effect  was 
f rontolytical  in  the  frontal  zone.  The  net  effect  of  both 
terms  yielded  high  frontogenetical  quantities  within  the 
frontal  zone  below  1000  feet  as  well  as  in  the  warm  air 
adjacent  to  the  front.  High  f rontolytical  values  were  found 
in  the  frontal  zone  above  1000  feet  where  the  "tilting"  term 
dominated  the  "confluence"  term. 

1.3.2  Coastal  Frontogenesis  Studies 

The  term  "coastal  frontogenesis"  was  first  used  by 
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Bosart (  1972  ) ,  however,  many  years  earlier  the  coastal  front 
concept  was  explored  in  relation  to  stratus  development  in 
the  southeast  United  States.  Carson (1950)  analyzed  three 
cases  of  "Gulf  Stream  f rontogenesis"  where  a  surface  warm 
front  developed  in  the  vicinity  of  the  Gulf  Stream's  central 
axis.  Carson  noted  the  synoptic  pattern  was  typically 
dominated  by  a  cold  core  high  pressure  wedge  and  that  the 
isobars  along  the  Atlantic  Coastal  Plain  were  oriented 
parallel  to  the  Gulf  Stream  axis.  West  of  the  warm  front, 
winds  were  northerly  with  much  cooler  temperatures  while 
wind  on  the  other  side  were  more  easterly.  Available 
soundings  indicated  that  air  over  the  land  at  1000-2000  feet 
was  10-15  K  warmer  than  at  the  surface. 

Carson's  studies  indicated  a  systematic  increase  of  the 
stratus  layer  tops.  The  stratus  layer  initiated  at  the  warm 
front  and  sloped  upward  toward  the  coast.  He  suggested  the 
stratus  layer  was  representative  of  the  slope  of  the  cold 
wedge.  Carson  hypothesized  that  air  flowing  over  the  warmer 
ocean  waters  was  lifted  and  carried  westward  by  the 
prevailing  winds  aloft.  The  strong  frontal  inversion 
trapped  the  lifted  air  making  it  possible  for  stratus  to 
form  and  develop  downward  as  a  result  of  evaporating 
precipitation . 

Bosart (1972)  was  credited  with  the  first  extensive 
mesoscale  analysis  of  a  coastal  front.  Special  attention 
was  given  to  the  24  December  1970  case  along  the  New  England 
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coast.  Analyses  revealed  temperature  gradients  as  great  as 
1  K/Km.  The  front  separated  a  stronger , easterly  and  more 
geostrophic  flow  over  the  ocean  from  a  lighter  northerly 
flow  over  land.  He  also  discovered  a  tendency  for  frontal 
stagnation  along  the  Boston-Providence  line  as  well  as  a 
precipitation  maximum  just  west  of  the  frontal  zone. 

Bosart  sites  geography,  friction,  land-sea  thermal 
contrast  and  latent  and  sensible  heating  as  the  causitive 
factors  of  frontogenesis.  Using  Petterssen ' s( 1956 ) 
development  equation,  he  argues  that  the  geography  of  New 
England  as  well  as  the  Carolinas  and  Texas  are  favorable  for 
this  type  of  frontogenesis  because  of  the  natural  tendency 
for  development  of  the  appropriate  vorticity.  In  addition, 
he  concludes  that  the  axis  of  maximum  precipitation  is 
consistent  with  the  quasi-geostrophic  omega  equation  since 
the  axis  is  within  the  region  of  maximum  lower  tropospheric 
warm  air  advection. 

As  an  extension  of  this  study,  Bosart (1975)  used  the 
same  case  to  analyze  the  effects  of  geostrophic  versus 
observed  deformation.*  Analyses  indicated  that  the  magnitude 
of  observed  deformation  surpassed  its  geostrophic 
counterpart.  Observed  deformation  fields  aliqned  well  with 
the  frontogenesis  region  while  the  geostrophic  fields 
displayed  a  preference  for  the  axis  of  the  anticyclonic 
wedge.  Also,  the  angle  of  the  isotherms  to  the  geostrophic 
deformation  axis  of  dilatation  far  exceeded  the  critical  45 


degrees  necessary  for  f rontogenesis.  He  concluded  that 
geostrophic  deformation  was  not  a  factor  in  coastal 
f rontogenesis. 

Using  57  coastal  front  cases,  Bosart  also  presented  a 
climatology  of  New  England  coastal  fronts.  As  in  the  case 
study,  the  precipitation  patterns  revealed  an  axis  of 
maximum  precipitation  just  west  of  the  preferred  coastal 
front  location,  while  the  thermal  patterns  identified 
December  as  the  month  of  greatest  land-sea  temperature 
contrast.  This  latter  fact  is  consistent  with  the  December 
maximum  of  coastal  front  occurrence. 

Further  investigation  of  the  precipitation  distribution 
associated  with  the  coastal  front  was  conducted  by  Marks  and 
Austin( 1979 ) .  Vertical  sounding  analyses  of  eight  coastal 
front  cases  consistently  identified  a  shallow  region  of 
convective  instability  at  600  mb  with  a  neutral,  well  mixed 
layer  above  it.  Analyses  of  radar  observations  indicated 
that  precipitation  was  oriented  in  bands  of  convective  cells 
which  apparently  originated  from  the  advancing  synoptic 
cyclone.  The  movement  of  these  cells  was  consistent  with 
the  winds  at  the  layer  of  convective  instability.  This 
pattern  of  convective  instability  associated  with  bands  of 
convective  cells  is  consistent  with  the  cyclone  model  of 
Harrold (  197  3 )  as  well  as  the  case  study  results  of 
Kreitzberg  and  Brown (1970). 

The  precipitation  analyses  of  Mark's  and  Austin 


reconfirmed  Bosart's  conclusion  concerning  higher 
precipitation  values  near  the  coastal  front  zone.  Since  the 
precipitation  bands  were  a  result  of  the  synoptic  cyclone, 
they  hypothesized  that  coastal  fronts  produce  a  local 
intensification  of  the  precipitation.  Calculations  of  water 
vapor  flux  divergence  revealed  a  flux  convergence  maximum 
over  the  frontal  zone.  Marks  and  Austin  conclude  that  the 
mechanism  for  precipitation  enhancement  is  the  creation  of 
low  cloud,  due  to  flux  convergence,  which  allows  for  the 
accretion  of  water  droplets  onto  the  hydrometeors  falling 
from  the  mesoscale  bands. 

Efforts  to  numerically  model  the  coastal  front 
phenomena  were  initiated  by  Ballentine ( 1980 ) .  He  used  a 
three-dimensional  primitive  equation  model  to  investigate 
the  effects  of  surface  friction,  latent  heating,  land-sea 
temperature  contrast,  synoptic  sclale  forcing  and  other 
physical  properties  upon  coastal  f rontogenesis.  The  "basic 
model",  which  included  all  the  physical  effects,  created  a 
typical  coastal  front  scenario.  Low-level  convergence 
developed  near  the  coastline  with  strong  vertical  motions  in 
the  warm  air  east  of  the  front. 

Ballentine  varied  the  model  structure  to  identify  the 
role  of  individual  physical  effects.  The  results  indicated 
latent  heat  release  and  surface  friction  were  not  primary 
factors  influencing  coastal  front  circulation,  however, 
land-sea  temperature  contrast  and  synoptic  scale  forcing 
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were  identified  as  significant  factors.  Reducing  the  sea 
surface  temperature  6  K  decreased  the  backing  of  the  winds 
over  the  land  resulting  in  lessened  convergence  and  vertical 
velocities  near  the  coast.  In  the  same  fashion,  removing 
the  effects  of  an  advancing  700  mb  short  wave  significantly 
reduced  the  circulation  within  the  coastal  front  zone. 

A  case  of  explosive  cyclogenesis  along  a  Carolina 
coastal  front  was  analyzed  by  Bosart ( 1981 ) . 

Three-dimensional  analyses  of  the  f rontogenesis  equation 
terms  produced  results  remarkably  similar  to  Sanders(  1955 ) . 
The  confluence  term  dominated  the  f rontogenet ical  effects. 
Vertical  cross  sections  of  winds  normal  to  the  front 
revealed  strong  lateral  shear  along  the  frontal  boundary 
with  significant  vertical  shear  above  the  front.  In 
addition,  a  zone  of  strong  vertical  motions  was  found  along 
the  warm  boundary  of  the  front.  Bosart  concluded  that 
cyclogenesis  initially  resulted  from  strong  low  tropospheric 
warm  air  advection  associated  with  the  coastal  front 
circulation.  The  front  then  acted  as  a  steering  mechanism 
for  the  cyclone  until  it  came  into  phase  with  an  intense 
short  wave  trough  moving  eastward  from  the  Ohio  Valley. 
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2.  Data  and  Methodology 

2.1  Surface  and  Upper  Air  Data 

Data  for  this  case  study  were  obtained  from  the 
National  Climatic  Center.  Hourly  observations  from  civilian 
and  military  reporting  stations,  for  the  period 
0000  GMT  21  -  0000  GMT  23  December  1983,  were  used  to 
construct  surface  analyses.  Wind  and  temperature 
observations  from  Coast  Guard  stations  and  offshore  buoys 
were  also  used  to  better  detect  the  development  and  onshore 
movement  of  the  coastal  front.  The  region  of  coverage 
included  the  Atlantic  coastal  states  from  Mew  Jersey  to 
Georgia  as  well  as  Tennesee,  Kentucky  and  West  Virginia. 

Cross  sections  were  constructed  from  radiosonde 
observations  within  the  same  region.  Precipitation  analyses 
were  constructed  from  six-hourly  data  on  airways 
observations  and  from  hourly  recording  station  networks. 

2.2  Objective  Analysis  Scheme  and  Grid  Contruction 
Generating  fields  of  kinematic  variables  as  well  as 

fields  of  gradients,  fluxes  and  frontogenesis  terms  required 
the  execution  of  numerous  finite-difference  calculations. 

To  simplify  the  process,  the  Barnes ( 1964 , 1973 )  objective 
analysis  scheme  was  used  to  create  gridded  fields  of  basic 
data.  This  particular  scheme  was  selected  for  its 
simplicity,  accuracy  and  cost  effectiveness. 

The  Barnes  scheme  uses  Gaussian  weight  functions  to 


generate  the  gndded  values.  Weights,  W  ,  are  assigned 

to  the  original  data,  f  (x,y),  based  solely  upon  the 

distance,  dn,  between  the  data  and  the  grid  point. 

Outside  of  a  specified  critical  radius,  rc ,  a  weight  of 

zero  is  assigned  to  the  data.  This  zero  weight  value  is  a 

local  modification  of  the  Barnes  scheme  which  reduces 

computation  time.  A  weight  parameter,  k  ,  is  externally 

choosen  to  determine  the  filter  response  function,  D  , 

o 

which  determines  what  percentage  of  an  original  wave's 
amplitude  is  retained  in  the  interpolated  field. 

With  one  pass  through  the  data,  the  filter  response 
function  can  be  shown  to  be 

Dq  =  exp (  (  kQ*2 jt/\  )  2  )  2.1 

where  is  a  given  wavelength.  The  values  at  the  grid 
points  are  then  calculated  using  Equation  2.2 
30  -  f„U,y)*Do  2.2 

Barnes(1973)  modified  the  original  technique  to  decrease 
computation  time.  The  revised  version  allows  for  a  second 


"correction  pass"  through  the  data  to  receive  the  desired 
response  for  a  qiven  wavelength.  This  is  acheived  by 
reducing  the  weight  parameter  on  the  correction  pass  to  a 


value  of  k^  where 


k.  =  y  *k 
1  o 


and  where  7  is  a  parameter  which  forces  the  interpolated 


field  to  converqe  toward  the  observed  field. 

The  computer  program  used  for  this  analysis  was 


designed  by  Dewey  Harms  of  this  department  as  a  part  of  the 

preliminary  studies  for  the  upcoming  GALE  project.  As  a 

part  of  his  project,  Harms(1985)  conducted  several  tests  on 

Barnes  scheme  results.  Analyses  of  data  created  from  a 

Midwest  cyclone  case  indicated  that  a  weight  parameter  of 
2  2 

250  Km  to  500  Km  yielded  the  minimum  Root  Mean 
Square  Error  values  for  most  meteorological  parameters. 

From  this.  Harms  determined  that  the  total  response 
function (a  measure  of  how  closely  the  interpolated  values 
match  the  observed  values)  for  a  wavelength  of  190  Km(two 
times  the  average  station  spacing)  was  0.936.  A  total 
response  function  value  of  1  indicates  an  exact  agreement. 

Using  this  information  and  other  guidelines  from  Harms, 
grid  regions  for  surface  and  upper  air  analyses  were 
developed.  Harms  suggested  that  grid  point  spacing  should 
range  from  one  half  to  one  times  the  average  spacing  of  data 
stations  within  the  grid.  The  nearest  neighbor  concept 
suggested  a  surface  station  spacing  over  the  region  of 
concern  was  about  85  Km.  The  value  of  60  Km  was  then 
selected  for  the  surface  grid  spacing.  Using  the  total 
response  function  calculated  by  Harms,  the  equations  of  the 
Barnes  scheme  were  solved  for  the  optimum  weight  parameter. 
A  value  of  400  km  was  derived  for  the  60  Km  spacing. 

An  upper-air  grid  was  developed  using  the  same 
procedures.  A  grid  spacing  of  150  Km  was  selected  based 
upon  a  280  Km  station  spacing.  An  upper  air  grid  weight 
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parameter  of  4600  Km  was  derived,  however,  the 

resulting  qridded  data  fields  were  clearly  too  "smooth".  A 

series  of  tests  using  variuos  weight  parameters  indicated 

2 

that  a  value  of  1000  Km  produced  reasonable  results. 

Figure  2.1  shows  the  domain  of  the  grids  used  in  this 
experiment.  The  surface  grid  extends  over  the  Atlantic  to 
include  the  available  buoy  data  and  remove  edge  effects  of 
the  scheme  away  from  the  coast.  Since  effectiveness  of  the 
Barnes  scheme  is  dependent  upon  the  selected  station 
spacing,  "bogus"  data  were  included  at  lxl  degree  locations 
over  the  ocean  where  the  data  void  would  adversely  affect 
the  results  in  the  generated  data  fields.  The  effect  of 
data  paucity  was  a  problem  only  in  the  southeast  sector  of 
the  upper-air  grid  where  one  "bogus"  sounding  was  added  to 
the  data.  The  "bogus"  data  were  developed  using  available 
surface  data  from  reporting  stations  and  off-shore  buoys, 
upper-air  charts  and  radar  film  loops.  In  most  cases,  a 
linear  interpolation  was  used  to  estimate  parameter  values. 
2.3  Frontogenesi s  Calculations 

Both  the  Eulerian  and  Lagrangian  perspectives  have 
distinct  advantages  when  quantifying  f rontoqenet ic 
processes.  The  Eulerian  method  focuses  upon  processes  which 
directly  affect  the  intensity  and  motion  following  the 
front.  Pettersen (  1940  )  defined  f rontogenesis  as  the  local 
derivative  of  the  gradient  of  some  conservative  property  S. 
Appendix  A  explains  how  equation  2.3  is  derived  from 
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Petterssen's  original  formula,  namely, 
d/dt  IVSl  =(  1/  IVSl  ) * { VS ) •  [V(  — Vh»  VS  )  +V( w«VS  )  +V ( dS/dt )  ]  2.3 

ABC 

The  formula  implies  that  f rontogenesis  is  a  function  of  the 
unit  gradient  vector  dotted  with  the  following:  the  gradient 
of  horizontal  S  advection  (term  A),  the  gradient  of  vertical 
S  advection  (term  B),  and  a  differential  diabatic  term  (term 
C).  For  our  purposes,  we  will  evaluate  only  term  A  and 
assume  that  terms  B  and  C  are  negligible.  This  reduces 
equation  2.3  to 

S/at  IVSl  =(  1/ IVSl  )*VS*V(V,  •  VS)  2.3.1 

h 

Figure  2.2  graphically  interprets  the  f rontogenesis 
processes  of  equation  2.3.1.  The  horizontal  gradient  vector 
is  directed  from  low  to  high  values  of  S(left  to  right).  To 
the  right  of  C,  the  horizontal  gradient  of  S  advection 
vector  is  directed  downstream,  while  to  the  left  of  C,  the 
horizontal  gradient  of  S  advection  vector  is  directed 
upstream(low  to  high  values).  From  equation  2.3, 
f rontogenesis  is  expected  to  occur  where  the  gradient  of  S 
advection  is  oriented  in  the  same  direction  as  the  S 
qradient  vector. 

The  Lagranqian  viewpoint  of  f rontogenesi s  focuses  on 
the  evolution  of  the  temperature  gradient  following  the 
motion  of  an  air  parcel.  Defining  frontogenesis  as  the 
total  derivative  of  the  gradient  of  the  conservative 
property  S,  Newton (1954)  derives  equation  2.4  as  the 


stration  of  effects  of  temperature  and  wind 
ds  upon  Eulerian  frontogenesis.  Wind  flow  is 


horizontal  components  of  the  frontoqenesis  equation  where 

the  equation  for  F  is  the  similar  to  F  . 

x  y 

Fy  =d/dt (dS/dy) = 

d/dy(dS/dt)-(dS/dx*du/dy+dS/ay*dv/ay+dS/dz*dw/dy)  2.4 
A  BCD 

Term  A  represents  the  effect  of  differential  diabatic 
heatinq  or  coolinq  and  is  difficult  to  qraphically  or 
numerically  evalute .  Terms  B,C  and  D  are  illustrated  in 
Figure  2.3.  Term  B  represents  the  effects  of  horizontal 
shear,  while  term  C  represents  the  effect  of  lateral 
confluence.  Term  D  illustrates  the  effect  of  the  horizontal 
gradient  of  vertical  velocity. 

Both  methods  of  assessing  frontoqenesis  have  distinct 
advantages.  As  previously  stated,  the  Laqrangian  method 
evaluates  the  f rontogenetical  effects  of  confluence,  shear, 
and  gradient  of  vertical  velocity  on  a  parcel  of  air  as  it 
flows  through  a  given  region.  On  the  other  hand,  the 
Eulerian  version  reveals  information  concerning  frontal 
movement.  For  example,  if  the  magnitude  of  the  qradient  is 
increasing  at  the  point  where  the  gradient  of  S  is  already  a 
maximum,  it  is  clear  that  there  is  frontal  intensification. 
Another  advantage  of  the  Eulerian  version  lies  in  the  need 
to  look  only  at  the  gradient  of  advection  in  relation  to  the 
temperature  iradient.  Both  of  these  quantities  are 
discernable  from  surface  or  upper-air  charts. 


2.4  Kinematic  Calculations 

Calculations  of  observed  divergence,  deformation  and 
relative  vorticity( about  the  vertical  axis)  are  performed 
usinq  the  standard  equations  of  Petterssen ( 1 940 ) .  Vertical 
velocities  are  obtained  usinq  the  kinematic  method.  The 
continuity  equation  in  X,Y,P  coordinates  is  expressed  as: 

du/dx  +  dv/dy  +  dw/dp  =0  2.5 

Integration  in  the  vertical  yields: 


(dw/dp)dp  =  Ve^V^  dp 


In  finite-difference  form,  the  vertical  velocity  at  a 

pressure  level  is  obtained  from  equation  2.7 

w  =w  -  V •  V,  *  ( p  -p.  )  2.7 

u  1  h  mj  *1 

where  w^  is  the  vertical  velocity  at 
the  upper  pressure  level 
w^  is  the  vertical  velocity  at 
the  lower  layer 

Pu  is  the  upper  pressure  level 

p^  is  the  lower  pressure  level 

V«V^  is  the  horizontal  divergence  of  the 

layer  mean  wind 


The  layers  used  to  calculate  layer  mean  divergences  were  the 
surface  to  975  mb  layer  and  successive  25  mb  layers  up  to 
800  mb  and  succesive  50  mb  layers  up  to  500  mb. 


3.  Analysis  of  the  Coastal  Front 

3.1  Synoptic  Overview 

Figure  3.1  illustrates  the  12-hourly  surface  and  500  mb 
charts  for  the  period  1200  GMT  21  December  throuqh  0000  GMT 
23  December  1983.  During  this  period  the  ridge  over  the 
eastern  U.S.  progresses  eastward  while  an  outbreak  of  arctic 
air  moves  southward  from  the  Northern  Plains  States.  An 
upper-air  cutoff-low  develops  and  advances  eastward  over  the 
northern  Great  Lakes  region. 

The  surface  pattern  shows  a  dome  of  high  pressure 
centered  over  New  England  and  eastern  Canada.  As  the  major 
axis  of  the  ridge  progresses  eastward,  an  anticyclonic  bulge 
develops  east  of  the  Appalachians.  By  0000  GMT  22  December, 
a  trough  develops  along  the  Mid-Atlantic  coastal  region  in 
association  with  a  well  defined  mesoscale  frontal  zone.  By 
1200  GMT,  the  trough  progresses  inland  and  a  weak  low  center 
evolves  near  the  North  Carolina-Virginia  border.  At 
0000  GMT  23  December,  the  synoptic  scale  cold  front 
associated  with  the  outbreak  of  cold  air  in  the  Midwest, 
progresses  into  the  western  Carolinas  and  begins  to 
overwhelm  the  mesoscale  features  associated  with  the  coastal 
front . 

3.2  Mesoscale  Analysis 

Figure  3.2  offers  a  detailed  mesoscale  perspective  of 
coastal  front  evolution.  Ship  data  was  rarely  available, 
however,  data  from  the  offshore  buoys  are  significant  in 


1200  GMT  21  Dec 


0000  GMT  22  Dec 


//( 


Figure  3.1. 
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Surface  and  500  mb  maps  for  21-23  December  1983. 
Standard  plotting  analysis  is  used.  Solid  lines 
represent  500  mb  heights,  in  decameters,  and 
surface  pressures,  in  millibars  in  excess  of 
1000.  500  mb  isotherms  are  in  °C.  National 

Weather  Service  products  were  used  to  develop  the 
analyses. 


identifyinq  and  tracking  the  coastal  front  evolution.  At 
1800  GMT  21  December  the  pressure  pattern  reveals  the 
developing  trough  over  the  Atlantic  and  the  increasingly 
ageostrophic  flow  over  the  land.  In  particular,  winds  near 
the  coast  are  nearly  geostrophic  while  at  the  axis  of  the 
anticyclonic  bulge  the  winds  are  nearly  perpendicular  to  the 
isobaric  pattern  and  parallel  to  the  isotherm  distribution. 
Since  the  flow  inland  results  in  weak  cold  air  advection 
while  the  coastal  flow  favors  significant  warm  air 
advection,  the  net  effect  of  continued  ageostrophy  would 
enhance  f rontogenesis  near  the  coast. 

At  0000  GMT  22  December  the  thermal  pattern  indicates 
an  increase  in  the  isotherm  packing  along  the  coastal  region 
of  the  Carolinas  and  Virginia.  Further  evidence  of  offshore 
trough  development  is  offered  as  the  buoy  data  reveals  a 
shift  to  southeasterly  winds.  Cape  Hatteras,  NC(HAT)  also 
indicates  a  gradual  shift  from  a  northeasterly  to  a 
southeasterly  flow  with  an  increase  in  the  sustained  and 
gust  wind  speeds.  The  anticyclonic  bulge  and  associated 
ageostrophic  flow  remains  firmly  entrenched  east  of  the 
mountains  maintaining  slight  cold  air  advection  down  the 
axis  of  the  ridge. 

The  0600  GMT  22  December  analysis  affords  the  most 
interesting  view  of  the  coastal  front.  In  particular,  winds 
all  along  the  coast  from  Charleston,  SC(CHS)  to  Atlantic 
City,  NJ(ACY)  have  shifted  to  the  southeast  with  a 


Figure  3.2.  (continued) 
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significant  increase  in  sustained  and  gust  wind  speeds.  The 
associated  trough  shifts  inland  significantly  and  appears  to 
sever  the  anticyclonic  bulge  over  central  North  Carolina  and 
western  Virginia.  The  total  wind  field  indicates  a  strong 
convergence  zone  along  coastal  South  Carolina  and  just 
inland  along  North  Carolina  and  Virginia.  Temperatures 
along  the  Carolina  coasts  are  increasing  at  modest  rates, 
however,  the  startling  temperature  contrast  lies  between  CHS 
and  the  Naval  Air  Station  at  Beaufort,  SC  A  temperature  of 
289K  (291K  the  previous  hour)  at  CHS  compared  with  a 
temperature  of  280K  at  Beaufort  converts  to  a  temperature 
gradient  of  approximately  lK/lOKm.  Considering  that  this 
gradient  is  oriented  somewhat  along  the  frontal  zone,  the 
true  cross  front  gradient  is  probably  much  greater. 

At  1200  GMT  22  December,  the  anticyclonic  bulge  which 
dominated  the  inland  flow  early  in  the  period  is 
non-existent.  The  coastal  front  moves  offshore  along  South 
Carolina  while  strong  southeasterly  winds  along  the  coast 
converge  with  northwesterly  winds  inland  to  increase  the 
temperature  gradient  along  the  eastern  third  of  North 
Carolina.  By  1800  GMT  22  December,  the  strong  temperature 
contrast  still  exists  across  Virginia  and  North  Carolina, 
however,  the  characteristic  convergent  flow  associated  with 
the  front  no  longer  exists.  A  weak  cyclone  center  located 
along  the  Carolina  coast  at  1200  GMT  progresses 
northeastward  and  deepens  through  the  1800  GMT  22  December 


time  period.  The  circulation  associated  with  this  cyclone 
results  in  a  westerly  flow  over  most  of  the  region.  The  net 
result  is  a  qradual  diffusion  of  the  temperature  gradient 
associated  with  the  coastal  front. 

Figure  3.3  illustrates  the  temperature  and  wind 
profiles  for  selected  stations  for  the  time  period  1800  GMT 
21  December  to  1800  GMT  22  December.  Analysis  of  these 
profiles  reveals  some  interesting  characteristics  of  the 
coastal  front.  At  CHS,  Pope  AFB ,  NC(POB)  and  Richmond, 
VA(RIC)  a  sharp  rise  in  temperature  over  a  one  to  three  hour 
period  is  preceded  by  a  shift  to  southeasterly  winds, 
however,  at  other  coastal  stations  like  Langley  AFB,  VA  and 
HAT  and  New  Bern,  NC(EWN)  the  shift  to  southeasterly  winds 
results  in  minor,  if  not,  negligible  temperature  rises.  The 
fact  that  inland  stations  experience  significant  temperature 
increases  with  the  shift  to  southeast  winds  while  some 
coastal  stations  experience  minimal  temperature  increases 
suggests  that  the  coastal  f rontogenesis  process  is  not 
always  limited  to  the  immediate  coastal  region.  Previous 
description  of  the  coastal  front  link  sensible  heating  of 
the  surface  layer  over  the  ocean  with  land-sea  differential 
friction  effects  to  create  conditions  favorable  for 
f rontogenesi s  at  the  coastline.  If  coastal  f rontogenesis 
always  occurs  in  this  manner,  we  should  expect  the 
temperature  profile  at  HAT  to  imitate  the  profile  at  CHS. 
Instead,  the  profile  of  POB  seems  to  resemble  the  CHS 
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pattern  the  closest.  At  CHS,  POB  and  even  RIC  the  equally 
sharp  temperature  variations  seem  to  indicate  that  the  front 
is  approximately  a  steady  state  phenomena  at  the  time  it 
passes  these  locations{ assuming  a  constant  rate  of 
movement).  In  contrast,  at  Seymour  Johnson  AFB,  NC(GSB)  and 
Myrtle  Beach,  SC(CRE)  moderate  temperature  increases 
associated  with  the  southeasterly  winds  are  followed  by 
smaller  temperature  increases  over  the  next  three  to  four 
hours.  At  these  same  stations  temperature  decreases,  in 
association  with  a  retreating  coastal  front,  are  much  more 
acute.  These  observations  suggest  that  atmospheric 
processes  favored  frontogenesis  along  the  coastline  east  of 
CHS,  however,  it  appears  that  frontogenesis  was  occurring 
over  the  land  in  eastern  North  Carolina. 

3.3  Kinematic  Analyses 

The  analysis  of  kinematic  variables  is  useful  for 
relating  the  characteristics  of  the  wind  field  to  the 
generation  and  maintenance  of  a  frontal  zone.  Figures  3.4 
through  3.7  display  fields  of  divergence,  relative 
vorticity,  deformation  and  sea  level  geostrophic  deformation 
at  six-hourly  time  periods  from  0000  GMT  to  1800  GMT  22 
December  1983.  Grid  spacing  for  all  finite  difference 
calculations  was  60  Km.  Values  over  the  open  ocean  must  be 
considered  suspect  since  they  are  mostly  generated  from  the 
"bogus"  data  mentioned  earlier. 

At  0000  GMT  the  divergence  and  relative  vorticity 


patterns  are  consistent  with  previous  coastal  front 
analyses.  A  broad  region  of  divergence  and  negative 
relative  vorticity  associated  with  the  anticyclonic  bulge, 
extends  along  the  Appalachian  chain.  The  positive  vorticity 
region  is  negligible  at  this  time,  however,  a  broad  region 
of  convergence  extends  along  the  coast.  The  convergence 
maximum  over  HAT  is  consistent  with  the  shift  to 
southeasterly  winds  at  that  location.  The  observed 
deformation  maximum  parallels  the  coast  while  the 
aeostrophic  deformation  is  located  well  inland.  These 
results  are  guite  similar  to  the  findings  of  Bosart ( 1975 ) . 
The  inland  maximum  of  geostrophic  deformation  aligns  well 
with  the  strong  anticyclonic  curvature  of  the  high  pressure 
bulge . 

By  0600  GMT  (Fig.  3.5)  the  relative  vorticity  and 
convergence  fields  shift  significantly  westward  as  does  the 
observed  deformation  field.  These  patterns  correspond  well 
with  the  surface  trough  which  moved  inland  at  this  time. 
Strong  southerly  winds  at  the  coastline  and  notherly  wind 
over  the  interior  translate  to  strong  convergence  values 
along  the  coastal  region.  Also,  the  vorticity  maximum 
coincides  with  the  nearly  closed  cyclonic  circulation 
between  Myrtle  Beach  and  CHS.  Considering  the  magnitudinal 
increases  of  these  fields  and  that  the  axis  of  dilatation 
associated  with  the  observed  deformation  field  is  oriented 
parallel  to  the  coast  and  the  isotherm  pattern,  it  is 
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reasonable  to  expect  an  increase  of  the  temperature  gradient 
over  the  coastal  region. 

The  1200  GMT  relative  vorticity  and  divergence  fields 
(Fig.  3.6)  undergo  a  significant  change.  The  major  axes  of 
the  positive  relative  vorticity  and  convergence  shift 
eastward  and  the  convergence  maximum  slides  northwards  to 
the  Chesapeake  Bay.  The  wind  pattern  confirms  the  shift  of 
the  convergence  maximum  as  northwest  winds  over  northeast 
Virginia  and  strong  southeasterly  winds  over  eastern 
Maryland  indicate  a  strong  convergent  flow.  The  geostrophic 
deformation  maximum  previously  located  along  the 
Appalachians  has  diminished  significantly.  The  diminishing 
geostrophic  deformation  pattern  correlates  well  with  the 
dissipation  of  the  anticyclonic  bulge.  The  observed 
deformation  field  begins  to  deteriorate  at  this  time  also. 
These  particular  conditions  provide  additional  evidence  that 
the  coastal  front  has  probably  reached  its  maximum  intensity 
and  has  made  its  farthest  inland  penetration. 

At  1800  GMT  (Fig.  3.7)  the  vorticity,  observed 
deformation  and  divergence  patterns  confirm  that  the 
Carolina  coastal  plains  region  is  no  longer  favorable  for 
f rontogenesis.  The  positive  vorticity  maximum  over  the 
ocean  is  especially  suspect  considering  that  the  cyclonic 
flow  is  located  farther  north.  The  convergence  maximum  over 
West  Virginia  is  undoubtedly  in  conjunction  with  the 
advancing  synoptic  scale  cold  front. 


Same  as  Figure  3.4  except  for  1800  GMT  22  December  1983 


3.4  Surface  Frontogenesi s  Calculations 

By  definition,  a  frontal  zone  is  a  region  of  relatively 
strong  temperature  gradient.  Changes  in  the  temperature 
gradient  pattern  should  reflect  changes  in  frontal  intensity 
and  location.  Fig.  3.8  illustrates  the  potential 
temperature  gradient  field  at  three-hourly  intervals  from 
0000  GMT  to  1500  GMT  22  December  1983.  The  patterns  display 
a  region  of  strong  temperature  gradient  along  the  coastline 
from  South  Carolina  to  Maryland.  Although  the  temperature 
gradient  maximum  shift  slightly  westward  between  0000  GMT 
and  0600  GMT  and  then  retreats  eastward  throughout  the 
remaining  periods,  no  significant  change  in  the  temperature 
gradient  magnitude  is  displayed.  These  results  give  little 
indication  of  significant  f rontogenesis  along  the  coastal 
region . 

3.4.1  Lagrangian  Frontogenetical  Analysis 

To  assess  the  surface  frontogenetical  features  from  a 
Laqrangian  perspective,  Equation  2.4  in  Section  2.2  is 
evaluated.  Centered  differencing  is  used  to  evaluate  all 
derivatives.  The  diabatic  term  and  the  "tilting”  term  are 
not  evaluated. 

Figure  3.9  displays  the  total  frontogenetical  field  at 
three-hourly  intervals  from  0000  GMT  to  1500  GMT  22 
December.  In  general,  the  shear  term  magnitudes  were  much 
less  than  that  of  the  confluence  term,  therefore  the 
frontogenetical  fields  for  the  individual  terms  are  omitted 


Figure  3.8.  Temperature  gradient  at  three-hourly  periods 
from  0000  to_1500(arf)  GMT  22  December  1983. 
Units  are  10  Km.  Contour  interval  is 
every  1  unit. 


Figure  3.9. 


Lagrangian ( total )  frontogenesis  for 
three-hourly  periods  from  0000  to  1500(a-f)  GMT 
22-??CembeI  l?83*  Units  are 
10  K  m  s  .  Positive  values 
indicate  frontogenesis.  Contour  intervals  are 
every  20  units. 


from  the  diagrams.  Units  of  f rontogenesis  can  be  multiplied 

g 

by  10  to  yield  units  of  approximately 
K  (10  Km)  (3Hr)  At  0000  GMT,  a  broad  region  of 

f rontogenesis  parallels  the  entire  coastal  region  with  the 
maximum  values  along  the  coast  of  the  Carolinas.  At  0300 
GMT  and  0600  GMT  the  frontogenet ical  field  indicates  a 
significant  increase  in  the  temperature  gradient  tendency 
over  eastern  South  Carolina  and  southeastern  North  Carolina. 
These  sharp  increases  reflect  the  broad  onshore  flow  in  this 
region  and  are  consistent  with  the  increasing  convergence 
values  along  the  coastal  regions. 

By  0900  GMT  the  frontogenesis  field  signals  a  rapidly 
increasing  temperature  gradient  tendency  over  eastern  North 
Carolina  while  along  the  South  Carolina  coast  the 
frontogenesis  values  have  decreased.  This  pattern  is 
consistent  with  the  eastward  movement  of  the  coastal  front 
away  from  CHS  and  the  maximum  inland  penetration  of  the 
front  in  North  Carolina.  At  1200  GMT  the  frontogenesis 
maximum  shifts  further  eastward  and  decreases  in  magnitude 
giving  indications  that  the  coastal  front  is  retreating.  At 
this  time,  the  winds  over  the  interior  are  predominantly 
from  the  north  while  the  coastal  winds  are  shifting  to  the 
south.  The  horizontal  shear  associated  with  this  flow 
pattern  is  consistent  with  maintaining  an  existing  frontal 
zone,  however,  it  is  not  as  conducive  to  increasing  frontal 
intensity.  This  may  account  for  the  decrease  in  the  overall 


f rontogenesis  magnitudes  and  the  fact  that  the  shear  and 
confluence  terms  are  now  about  equal  in  magnitude.  The 
f rontogenesis  field  loses  further  definition  at  1500  GMT  as 
a  southwesterly  flow  begins  to  dominate  the  entire  region. 
3.4.2  Eulerian  Frontogenetical  Analysis 

The  Eulerian  frontogenesis  calculations  for  this 
analysis  are  based  upon  Equation  2.3  in  Section  2.2.  The 
frontogenesis  values  are  a  function  of  the  unit  gradient 
vector  dotted  with  the  gradient  of  advection.  In  essence, 
the  results  of  these  calculations  represent  the  effects  of 
the  thermal  advection  field,  derived  by  the  Barnes  analysed 
data,  upon  the  local  change  of  the  temperature  gradient. 

g 

Again,  the  values  can  be  multiplied  by  10  to  obtain 
units  of  approximately  K  (lOKm)  (3Hr)-1. 

Figure  3.10  shows  the  Eulerian  frontogenesis  patterns 
at  three-hourly  time  periods  from  0000  GMT  to  1500  GMT  22 
December.  At  0000  GMT  and  0300  GMT  a  weak  region  of 
frontogenesis  over  the  coastal  region  of  South  Carolina  is 
co-located  with  the  temperature  gradient  maximum,  in  Fig. 
3.8,  along  the  coastal  region  of  the  Carolinas.  As 
described  earlier  this  pattern  suggests  an  increase  in  the 
frontal  zone  temperature  gradient  without  significant 
frontal  movement  over  South  Carolina  and  slight  westward 
movement  over  North  Carolina.  The  configuration  of  the 
frontogenesis  field  changes  at  0600  GMT  as  the  maximum 
shifts  westward.  Considering  the  previous  location  of  the 


Fiqure  3.10.  Eulerian  f rontogenesis  for  three-hourly 

periods  from  0000  to  1500 (a-f)  GMT  22  December 
1983.  Units  and  contour  interval  are  as  in 
Figure  3.9.  Positive  values  indicate 
f rontogenesi s . 


temperature  gradient  maximum,  the  pattern  predicts  a 
westward  movement  of  the  frontal  zone.  The  magnitude  and 
location  changes  in  the  f rontogenesis  field  physically 
represent  the  increasing  coastal  warm  air  advection  clashing 
with  the  cold  air  advection  developing  over  western  South 
Carolina  and  resulting  in  a  larger  gradient  of  temperature 
advection . 

At  0900  GMT  the  f rontogenesi s  maximum  shifts  to  eastern 
North  Carolina  as  the  cold  air  advection  sweeps  eastward 
across  South  Carolina  and  into  central  North  Carolina.  This 
illustration  verifies  the  offshore  frontal  movement  along 
the  South  Carolina  coast  and  frontal  intensification  with 
slight  eastward  movement  over  North  Carolina.  The  1200  GMT 
and  1500  GMT  pattern  indicate  further  eastward  movement  and 
decrease  in  magnitude  of  f rontogenesis.  This  is  a  direct 
response  to  a  westerly  flow  which  develops  over  the  region, 
resulting  in  the  extension  of  cold  air  advection  all  the  way 
to  the  coastline  of  North  Carolina. 

Figure  3.11  represents  the  actual  local  rate  of  change 
of  the  temperature  gradient,  as  calculated  from  the  Barnes 
analyzed  fields,  for  0900  and  1200  GMT.  These  values  are 
obtained  by  subtracting  the  temperature  gradient  values  at 
four-hour  intervals  (two  hours  on  either  side  of  the 
specified  time  period).  Comparing  these  fields  with  the 
Eulerian  f rontogenesis  fields  in  Figure  3.10  it  is  clear 
that  the  Eulerian  calculations  are  predicting  too  large  a 


Fiqure  3.11.  Actual  local  change  of  the  temperature  gradient 
for  0900  GMT(a)  and  1200  GMT(b)  22  December. 
Units  are  as  in  Figure  3.9.  Contour  interval  is 
every  40  units.  Positive  values  indicate  an 
*  increase  in  the  temperature  gradient. 


temperature  gradient  increase  in  the  frontal  zone.  Another 
interesting  feature  is  that  both  the  Lagrangian  and  Eulerian 
equations  predict  frontogenesis  yet  the  temperature  gradient 
patterns ( Fig .  3.8)  do  not  show  any  significant  change.  The 
time  limit  over  which  these  values  are  calculated  may  be  the 
cause  of  this  discrepancy.  The  frontogenesis  calculations 
are  instantaneous  values,  however,  the  actual  change  in  the 
temperature  gradient  is  a  result  of  a  four-hour  time 
differential  calculation.  Considering  that  the  front  moved 
inland  and  then  retreated  during  this  time,  a  major  shift  in 
the  temperature  gradient  maximum  probably  would  not  appear. 
3.5  Vertical  Structure  of  the  Coastal  Front 

Figure  3.12  illustrates  vertical  soundings  for  several 
locations  in  the  region  of  the  coastal  front.  The  soundings 
of  stations  on  the  cold  air  side  of  the  front  are 
characterized  by  a  strong  temperature  inversion  below 
900  mb.  For  example,  the  1200  GMT  22  December  sounding  at 
Charleston  reveals  a  6.5  K  temperature  increase  between  the 
989  mb  and  974  mb  levels.  Calculations  using  the 
hypsometric  relationship  convert  this  to  approximately  a 
1  K/20  meter  vertical  temperature  gradient.  The  cold  wedge 
clearly  creates  a  highly  stable  lower  boundary  layer  which 
partially  accounts  for  the  persistent  nature  of  this 
phenomena . 

The  potential  temperature  cross  section  in  Figure  3.13 
provides  additional  information  about  the  three-dimensional 
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Figure  3.12.  Stuve  diagrams  of  temperature  and  dew  point 
at  Cape  Hatteras  and  Charleston  for21200  GMT 
22  December  1983.  Pressure  in  1(j  mb  and 
temperature  in  °C. 
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framework  of  this  coastal  front.  Although  the  overall 
pattern  suggests  generally  warmer  temperatures  along  the 
coastal  region  in  contrast  to  the  temperatures  west  of  the 
Appalachians,  the  mesoscale  structure  of  the  coastal  front 
is  the  most  prominent  feature  of  the  analysis.  The  thermal 
structure  of  the  front  is  characterized  by  the  intense 
horizontal  temperature  gradient  midway  between  Cape  Hatteras 
and  Greensboro.  This  intense  temperature  gradient  extends 
vertically  up  to  approximately  the  925  mb  level.  The 
maximum  depth  of  the  cold  air  appears  to  exist  at  or  near 
Greensboro.  This  is  consistent  with  the  previous  location 
of  the  anticyclonic  bulge  in  the  surface  pressure  pattern 
since  the  anticyclonic  bulge  is  a  physical  representation  of 
the  depth  of  the  cold  air. 

Eguation  2.4  (Section  2.2)  is  used  to  evaluate  the 
three-dimensional  f rontogenetical  patterns  associated  with 
the  coastal  front.  All  terms  except  the  diabatic  term  are 
estimated.  Gridded  fields  for  each  term  were  derived  using 
the  grid  shown  in  Fig.  2.1.  Values  were  calculated  at  the 
same  levels  specified  in  Section  2.4.  Values  were  then 
transposed  to  the  cross  sections.  Frontogenetical  values 
were  computed  using  0000  GMT  22  December  data,  however,  the 
frontogenetical  patterns  revealed  little  in  comparison  with 
the  1200  GMT  data. 

Figure  3.14a  -  3.14d  illustrate  the  frontogenetical 
fields  for  the  HAT  to  Nashville,  TN(BNA)  cross  section. 


BNA  GSO  HAT  . 

Figure  3.14.  Contribution  to  Lagrangian  frontogenesis  from 

confluence (a ) ,  shear(b),  and  tilting(c)  terms  and 
sum(d)  of  a,  b  and  c  at  1200  GMT  22  December 
1983.  Dashed  lines  are  276  K  and  288  K  potential 
temperature  lines. 


Horizontal  confluence  (Figure  3.14a)  and  horizontal  shear 
(Figure  3.14b)  are  frontogenetical  within  the  frontal  zone 
with  maximum  confluence  values  just  above  the  950  mb  level 
and  the  maximum  shear  contribution  just  below  950  mb.  The 
horizontal  gradient  of  vertical  velocity  ("tilting")  term  is 
predominately  f rontolytical  within  the  frontal  zone  with  the 
maximum  frontolytic  contribution  located  west  of  the 
f rontogenesis  maximums  of  the  shear  and  confluence  terms. 
Although  the  magnitudes  are  not  the  same,  these  results  do 
resemble  the  frontogenesis  patterns  of  the  coastal  front 
study  by  Bosart(1981)  and  of  the  intense  surface  cold  front 
study  by  Sanders(  1955  ) . 

The  wind  components  perpendicular  and  parallel  to  the 
HAT  to  BNA  cross  section  are  shown  in  Figure  3.15.  The 
patterns  display  a  strong  confluence  and  shear, 
respectively,  at  the  950  mb  level.  Within  the  coastal  front 
zone,  the  winds  are  predominately  from  the  north  and  west, 
while  the  winds  outside  the  frontal  zone  are  mostly  from  the 
south  and  east.  Particularly  striking  is  the  manner  in 
which  the  winds  perpendicular  to  the  cross  section  slope 
with  and  give  additional  definition  to  the  top  of  the 
frontal  zone.  The  frontal  zone  clearly  separates  the  colder 
northerly  flow  from  the  warmer  southerly  flow. 

The  strongly  convergent  flow  at  the  eastern  edge  of  the 
frontal  boundary  suggests  strong  vertical  motions  within  the 
.  Vertical  velocities  were  derived  using  the 
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kinematic  method  outlined  in  Section  2.3.  Figures  3.16  and 
3.17  illustrate  the  vertical  velocity  patterns  for  0000  GMT 
and  1200  GMT  22  December  along  the  HAT  to  BNA  cross  section 
as  well  as  the  Wallops  Island,  VA(WAL)  to  Dayton,  OH(DAY) 
cross  section.  Both  cross  sections  indicate  negligible 
vertical  velocities  along  the  coastal  region  at  0000  GMT, 
however,  the  1200  GMT  patterns  display  a  markedly  different 
configuration.  Ascent  of  about  2-3  x  10  ^  mb  s  * 
extends  into  the  upper  levels  of  the  coastal  front  zone. 

The  ascending  motion  at  the  lower  levels  is  apparently  an 
extension  of  the  ascending  motion  located  at  higher  levels. 
This  suggests  that  the  vertical  motions  within  the  zone  is 
not  solely  caused  by  the  convergence  within  the  frontal 
zone,  but,  may  be  characteristic  of  the  broader  scale  low 
and  mid-tropospheric  flow  pattern.  These  results  are 
somewhat  contradictory  to  the  results  of  Bosart(1981)  in 
which  a  separate  and  distince  region  of  strong  ascent,  in 
excess  of  12  x  10  ^  mb  s  ^ ,  within  the  frontal  zone. 

Although  the  ascending  motion  is  not  as  distinct  as  in  the 
case  investigated  by  Bosart,  the  general  circulation  within 
this  coastal  front  is  characteristic  of  the  thermally  direct 
circulation  discussed  by  Bosart.  The  air  in  the  warm  sector 
rises  as  it  passes  above  the  cold  air  and  the  cold  air  at 
the  surface  flows  back  toward  the  warm  sector. 

3.6  Precipitation  Analysis 

Fig.  3.18  illustrates  the  six-hourly  precipitation 


Figure  3.17.  Vertical  motions  along  the  Wallops  Island  to 
Dayton  cross  section  for  0000  GMT ( a )  and 

1200  GMT(b)  22  December.  Units  are  same  as 
Figure  3.16. 
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totals  for  designated  periods.  An  orderly  progression  of 
the  precipitation  maximum  from  western  North  Carolina  to 
southwestern  New  Jersey  is  evident.  The  maximum  over 
western  North  Carolina  probably  reflects  some  orographic 
effects.  Northeast  and  easterly  winds  through  0600  GMT 
undoubtably  resulted  in  some  lifting  up  the  mountain  slopes 
and  higher  precipitation  totals.  This  assumption  is  further 
substantiated  by  the  1200  GMT  totals  which  indicate  little 
precipitation  over  the  same  area  when  northwesterly  winds 
dominated . 

The  advancement  of  the  large  scale  precipitation 
maximum  is  consistent  with  the  vertical  velocity  patterns 
illustrated  in  Figs.  3.16  and  3.17.  The  strongest  vertical 
velocities  are  between  BNA  and  Greensboro,  NC(GSO)  at  0000 
GMT  22  December,  however,  the  vertical  velocity  maximum 
shifts  northeastward  to  the  vicinity  of  Washington,  DC(IAD) 
by  1200  GMT.  The  high  correlation  between  the  precipitation 
and  vertical  velocity  maximums  instills  confidence  that  the 
vertical  velocity  computations  are  representative  of  the 
atmospheric  motion. 

Bosart( 1972, 1975, 1981 )  and  Marks  and  Austin(1979) 
identify  the  cold  air  side  of  the  coastal  front  zone  as 
favorable  for  precipitation  enhancement.  Figure  3.18  seems 
to  indicate  similar  results.  In  Fig.  3.18(a),  a  maximum 
extends  through  central  South  Carolina  while  in  Fig. 

3.18(b),  a  maximum  is  evident  from  eastern  Virginia 
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southward  through  east  central  North  Carolina.  Both  of 
these  maxima  are  located  near  the  edge  of  the  farthest 
inland  penetration  of  the  coastal  front. 

Radar  observations  were  also  reviewed  to  study  the 
characteristics  of  the  precipitation  field.  Radar  films  and 
recorded  observations  from  CHS,  Wilmington,  NC(WIL),  HAT  and 
Patuxent  River,  MD(NHK)  were  analyzed.  In  an  attempt  to 
identify  precipitation  enhancement,  individual  cells  and 
band  movements  were  analyzed  in  conjunction  with  hourly 
precipitation  reports.  Inconsistencies  in  reporting 
accuracies  and  location  of  precipitation  in  relation  to 
reporting  stations  prevented  the  development  of  any 
conclusions  concerning  the  enhancement  of  precipitation  as  a 
cell  or  band  passed  through  the  coastal  front  zone. 

At  0000  GMT  22  December  1983,  the  CHS  radar  indicated  a 
broad  region  of  light  precipitation.  The  echo  pattern  was 
characteristic  of  stratiform  precipitation  with  few 
intensities  above  3  mm  hr  Between  0440  GMT  and  0520 

GMT,  a  broken  line  of  echoes  developed  into  a  solid  line  of 
showers  with  embedded  convective  activity.  The  line 
extended  100  Km  southward  from  a  point  20  Km  east  of  CHS. 

Fig.  3.19  illustrates  how  this  line  progressed  eastward  and 
embedded  convective  cells  developed.  The  band  was 
apparently  frontal  activity  associated  with  the  eastward 
movement  of  the  convergence  zone  associated  with  the  coastal 


front 


3.19.  Simulation  of  radar  PPI  scope  at  Charleston 
at  0530  GMT ( top )  and  0630  GMT(bottom) 

22  December.  VIP  levels  are  illustrated  as 
follows:  1( hatched),  2( cross-hatched ) ,  3 (dotted) 
4 ( black ) . 


Analysis  of  the  WIL  and  HAT  radars  indicates  small 
cells  of  convective  activity.  The  pattern  was  somewhat 
random  and  had  no  correlation  with  the  inland  movement  of 
the  coastal  front.  However,  a  small  line  of  showers  did 
develop  around  0700  GMT,  southwest  of  WIL,  and  progressed 
eastward.  This  was  probably  an  extension  of  the  line  which 
developed  east  of  CHS  a  few  hours  earlier. 

The  most  interesting  radar  observations  were  at  NHK. 
Several  distinct  bands  developed  over  eastern  Virginia  and 
progressed  northeastward  at  about  20  m  s-1.  The  bands 
often  expanded  as  they  moved  northeastward  and  frequently 
merged  with  other  areas  of  precipitation  in  the  northern  and 
eastern  region  of  radar  coverage.  Figure  3.20  illustrates  a 
few  of  these  bands.  Associated  rainfall  rates  were  on  the 
order  of  10-20  mm  hr  1 ,  which  is  consistent  with  the 
larger  vertical  velocities  and  observed  rainfall  amounts  in 
this  region. 

A  conection  between  mesoscale  rainbands  like  those 
identified  on  the  NHK  radar  and  convective  instability  was 
established  by  Kreitzberg  and  Brown (1970).  Marks  and  Austin 
used  those  concepts  to  link  rain  bands  with  coastal  front 
precipitation.  Fiqure  3.21  presents  vertical  profiles  of 
equivalent  potential  temperature  for  HAT,  GSO,  CHS,  and  IAD 
respectively.  The  GSO  and  CHS  temperature  profiles  are 
quite  similar  to  the  profiles  presented  by  Marks  and  Austin. 
An  inversion  layer  is  capped  by  a  deep  layer  of  warmer,  more 
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21.  Vertical  profiles  of  equivalent  potential 

temperature  at  0000  GMT  (solid)  and  1200  GMT 
(dashed)  22  December  at  Washington,  D.C.(a),  Cape 
Hatteras(b),  Greensboro ( c )  and  Charleston (d ) . 


stable  air.  Above  the  warmer  air,  a  convectively  unstable 
layer  of  about  100  mb  depth  is  evident.  The  HAT  profile 
displays  different  characteristics  since  this  station  is 
within  the  warm  air.  Warm  air  extends  up  to  the  650  mb 
level  where  a  convectively  unstable  layer  begins  and  extends 
to  the  580  mb  level.  The  convectively  unstable  layers  at 
IAD  are  much  more  shallow  and  are  of  lesser  magnitude. 

At  1200  GMT,  the  unstable  layers  have  mostly  decreased 
in  depth  and  magnitude  although  distinct  unstable  layers 
clearly  exist  at  HAT,  GSO  and  CHS.  On  the  contrary,  the 
unstable  layer  at  IAD  was  very  shallow( about  20  mb).  In 
studying  these  profiles,  it  is  evident  that  the  HAT,  GSO  and 
CHS  unstable  layers  are  deeper  and  apparently  more  unstable 
than  the  layer  at  IAD,  yet  the  precipitation  totals  were 
significantly  higher  at  the  latter  station.  These  profiles 
suggest  that  greater  amounts  of  precipitation  might  be 
expected  at  the  more  southern  stations;  however,  upward 
vertical  motion  must  exist  at  the  level  of  instability  in 
order  for  it  to  have  any  effect  upon  the  vertical 
velocities.  Figures  3.16(b)  and  3.17(b)  indicate  that 
little  upward  vertical  velocity  exists  at  the  instability 
level  at  three  southern  stations  while  at  IAD  strong  rising 
motion  exists  at  the  level  of  instability.  It  is  possible 
that  the  upward  vertical  motion  at  IAD  leads  to  the 
realization  of  the  covective  instability  thus  creating  a 
more  shallow  layer  of  instability. 
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Fig.  3.22  illustrates  the  wind  components  parallel  and 
perpendicular  to  the  movement  of  the  precipitation  fields  at 
HAT  and  I AD .  The  v-component  (perpendicular  to  the  band 
movement)  has  a  value  of  zero  at  the  steering  level  of  the 
precipitation.  At  both  stations  the  steering  level  is 
associated  with  the  base  of  the  convectively  unstable  layer. 
These  results  are  consistent  with  the  previous  findings  of 
Marks  and  Austin  and  suggests  a  relationship  between  the 
convective  instability  and  the  precipitation  fields. 


Wind  Speed 

Figure  3.22.  Wind  components  perpendicular (dashed )  and 

parallel ( solid )  to  precipitation  band  movement 
for  Cape  Hatteras  (filled  circle)  and  Washington 
D.C.( square)  at  1200  GMT  22  December.  Wind 
speeds  in  units  of  m  s~  . 


4. 


Geostrophic  Winds  and  the  Cold  Wedge 
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4.1  Alternative  Geostrophic  Wind  Computation  Methods 

The  horizontal  pressure  pattern  is  of  the  utmost 
importance  in  producing  horizontal  atmospheric  motions, 
however,  in  mountainous  terrain,  sea  level  pressure  patterns 
may  not  accurately  depict  the  correct  pressure  gradients. 

The  errors  are  largely  a  result  of  the  twelve-hour  mean 
temperature  factor  used  to  reduce  pressure  values  to  sea 
level  and  the  empirical  terms  added  to  pressures  above  1000 
ft.  Saucier ( 1955 )  indicates  that  errors  may  be  as  great  as 
lOmb  at  high  elevations  while  neighboring  stations  above  and 
below  1000  ft.  may  have  pressure  differences  as  much  as 
0 . 5mb . 

These  pressure  errors  often  result  in  erroneous  surface 
geostrophic  wind  calculations.  To  combat  this  problem, 
Sangster ( 1960 )  developed  a  new  method  to  compute  the 
pressure  qradient  force  on  a  non-isobaric  surface.  The 
method  is  based  upon  a  reference  surface  which  is  a  highly 
smoothed  representation  of  the  earth's  terrain.  Using  this 
surface  requires  the  reduct  ion ( both  upward  and  downward)  of 
surface  pressures  to  the  reference  level.  Sangster  claimed 
this  method  reduced  the  errors  since  most  pressure 
reductions  were  minimal  when  compared  to  the  pressure 
reductions  to  sea  level  from  high  elevations. 


Sanqster  expressed  the  pressure  qradient  force  as  an 
additive  function  of  a  geostrophic  stream  function  and  a 


geostrophic  potential  function.  His  test  results  indicated 
that  the  smoothed  terrain  geostrophic  winds  in  mountainous 
regions  were  often  more  accurate  than  the  geostrophic  winds 
assessed  from  the  sea  level  pressure  pattern.  His  method 
was  most  effective  in  synoptic  situations  when  cold  air 
bordered  along  a  mountain  range.  In  these  scenarios, 
extreme  sea  level  isobar  concentration  suggested  strong 
geostrophic  winds  yet  the  actual  winds  were  relatively 
light . 

Richwein ( 1980 )  examined  the  ineffectiveness  of  sea 
level  isobars  to  reflect  the  wind  flow  in  a  cold  wedge  east 
of  the  Appalachians.  He  employed  a  technique  for 
calculating  surface  geostrophic  winds  based  upon  the  thermal 
wind  relationship,  where  the  thermal  wind  is  the  geostrophic 
wind  at  a  lower  level  subtracted  from  the  geostrophic  wind 
at  an  upper  level.  Assuming  the  850  mb  level  to  be  a 
suitable  upper  geostophic  level,  Richwein  subtracted 
calculated  thermal  winds  from  the  850mb  winds  to  obtain 
surface  geostrophic  winds.  The  results  are  remarkably 
similar  to  the  actual  surface  winds  in  a  cold  wedge. 

Richwein  then  related  this  theory  to  coastal 
f rontogenesis.  Prior  to  f rontogenesis,  the  horizontal 
temperature  gradient  is  relatively  small  resulting  in  a 
small  thermal  wind.  Adding  the  thermal  wind  to  the  surface 
wind  produced  a  light  northeast  850  mb  wind.  This  result  is 
consistent  with  the  early  period  of  a  cold  wedge.  At  the 


mature  stage  of  f rontogenesis,  the  horizontal  temperature 
gradient  has  increased  significantly  resulting  in  a  strong 
thermal  wind.  Vectorially  adding  this  thermal  wind  to  the 
surface  wind  results  in  a  southerly  850  mb  wind.  This 
southerly  850  mb  wind  is  consistent  with  the  veering  of 
850  mb  winds  during  the  onset  of  frontogenesis. 

4.2  A  Variation  of  the  Sangster  Model 

Geostrophic  wind  computations  within  the  cold  wedge  are 
generally  poor  estimators  of  the  actual  wind  flow. 
Bosart(1975)  attributes  the  ageostrophy  to  differential 
friction  effects.  Richwein  and  Sangster  argue  that  sea 
level  pressure  patterns  do  not  accurately  depict  the 
horizontal  pressure  gradient  in  mountainous  terrain.  To 
further  investigate  this  problem,  modifications  are  made  to 
the  Sangster  model  to  test  the  hypothesis  that  pressure 
reductions  are  partially  responsible  for  the  apparent 
ageostrophy  of  the  cold  wedge  wind  pattern. 

Equation  4.1  expresses  the  geostrophic  wind  normal  to 
the  isobar  gradient. 

Vg=-(1/Pf )*dP/dn  4.1 

where  p=Density 

f=coriolis  parameter 
dP/^n=horizontal  pressure  gradient 
As  mentioned  earlier,  the  use  of  sea  level  pressures,  in 
equation  4.1,  can  lead  to  errors  in  mountainous  terrain. 
Sangster  avoided  this  problem  by  creating  a  fictitious, 


smoothed  surface.  Another  option  is  to  change  the  reference 
surface  from  mean  sea  level  to  a  terrain  level  reference 
surface.  By  usinq  the  earth's  terrain  as  a  reference  level, 
pressure  reductions  are  totally  avoided.  This  reference 
level  is  developed  by  calling  the  terrain  surface  <T  defined 
by : 

<T  =  P/P 

s 

Here,  P=pressure  at  a  given  level 

Ps=station  pressure 

One  sees  that  er  =  l  at  the  earth's  surface  and  is  the 

reference  surface  aforementioned. 

The  horizontal  pressure  gradient  on  a  cr-surface  is 

mathematically  derived  through  the  use  of  the  transformation 

equation  expressed  in  Equation  4.2. 

V  (Q)=V  (Q)+  0/  z*(V  (z))  4.2 

t  z  r 

Using  the  new  reference  surface,  where  r=  and  Q=P,  changes 
Equation  4.2  into  Equation  4.3. 

Va(P)=Vz(P)-  p/  z*(Vff(z))  4.3 

Assuming  a  hydrostatic  balance.  Equation  4.3  becomes 
(l/p)Vz(P)  =  (l/p)V<7  (P)+gVa(z)  4.4 

where  g=acceleration  of  gravity. 

The  left  hand  side  of  4.4  is  the  horizontal  pressure 
gradient  force  on  a  constant  height  surface  and  is  a 
function  of  the  gradients  of  station  pressure  and  elevation 
on  a  CT-surf ace. 

In  finite  difference  form,  the  horizontal  pressure 


4.5 


(1/p) ( dP/dn )  =( 1/p ) (dP/dn )  +g ( d z/dn ) 

Z  (To 

Substituting  Equation  4.5  into  Equation  4.1  yields  the 
geostrophic  wind  on  a  o-surface. 

V  ((r)»-(l/f)  (  (  1/p)  P/  n+g  z/  n)^  4.6 

4.3  Test  Results  of  the  Modified  Sangster  Method 

The  Barnes  analysis  scheme  is  used  to  create  gridded 
values  of  station  pressure  and  elevation.  Since  the 
pressure  gradient  is  a  function  of  both  variables, 
consistent  analyses  demand  that  elevations  are  used  for  only 
stations  which  have  available  station  pressures.  This  is  a 
minor  disadvantage  to  this  method,  since  a  new  gridded 
elevation  field  is  computed  for  every  hourly  data  set. 
Considering  the  minor  cost  of  the  additional  computations, 
this  method  is  still  much  preferred  to  Sangster' s  more 
costly  method  of  arbitrarily  reducing  pressures  to  a 
reference  surface. 

Figure  4.1  represents  the  Barnes  gridded  wind  vectors 
at  12Z  21  December  for  surface  winds,  sea  level  geostrophic 
winds  and  terrain  level  qeostrophic  winds  respectively. 

Both  geostrophic  wind  fields  display  similar  patterns  and 
are  in  sharp  contrast  to  the  actual  wind  field.  In 
particular,  two  discrepancies  are  readily  visible.  First, 
the  geostrophic  patterns  indicate  easterly  flow  over  the 
coastal  regions  while  northeast  flow  dominates  the  actual 
wind  pattern.  Second,  the  geostrophic  winds  quickly  shift 
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Figure  4.1.  Barnes  analyzed  wind  fields  for  actual  winds(a), 
sea  level  geostrophic  winds(b)  and  terrain  level 
geostrophic  winds(c)  at  1200  GMT  21  December. 
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to  the  south  over  the  more  elevated  terrain  while  the  actual 
wind  field  gradually  shifts  to  the  east  and  decreases  in 
magnitude.  These  discrepancies  are  consistent  with  the  sea 
level  isobar  pattern  which  is  oriented  east-west  over  the 
coastal  region  and  is  tightly  packed  in  a  north-south 
orientation  against  the  mountains. 

A  more  rigorous  examination  of  the  geostrophic  fields 
reveal  some  minor  variations  between  the  patterns.  The 
terrain  following  geostrophic  winds  are  slower  to  shift  to  a 
southerly  component  in  the  vicinity  of  the  anticyclonic 
axis.  This  variation  may  be  a  result  of  the  gradient  of 
elevation  term.  Over  the  coastal  regions,  the  gradient  of 
elevation  is  small  and  contributes  little  to  ( <r) , 
however,  along  the  mountain  chain,  the  elevation  gradient 
term  is  large  and  therefore  has  a  greater  impact  upon  the 
resultant  winds. 

To  test  the  accuracy  of  both  geostrophic  winds  in 
relation  to  the  actual  wind  fields.  Standard  Errors(SE)  are 
calculated  over  the  entire  region  from  the  Appalachian 
mountains  eastward.  SE's  are  also  calculated  for  just  the 
mountainous  region  to  determine  if  the  V  (o-) 
calculations  are  more  accurate  over  higher  sloped  terrain  as 
is  indicated  by  visual  inspection.  SE  values  were  obtained 
by  summing  the  difference  of  the  specific  geostrophic  wind 
value  from  the  Barnes  gridded  actual  wind  value.  This 
particular  method  requires  the  assumption  that  the  Barnes 


gridded  wind  field  represent  the  true  wind  field  over  the 
grid  region.  The  results  are  tabulated  in  Table  4.1  and 

4.2. 

Standard  error  values  indicate  that  V  ( o-) 

g 

calculations  over  the  highly  sloped  terrain  are  slightly 
more  accurate  than  over  the  region  as  a  whole.  The 
u-component  (north-south)  errors  for  terrain  level  winds  are 
as  much  as  5.5  m/s  less  than  the  geostrophic  winds  over  the 
sloped  terrain.  Over  the  region  as  a  whole,  the  u-component 
errors  only  differed  by  1-3  m/s.  Of  much  greater 
significance  are  the  overall  magnitudes  of  the  SE ' s  which 
range  from  6-15  m/s.  Since  winds  in  the  cold  wedge  rarely 
exceed  7  m/s,  it  is  clear  that  neither  method  is 
satisfactory  for  predicting  the  observed  wind  pattern. 

To  insure  that  density  calculations  were  not 
influencing  the  wind  computation,  density  values  were 
obtained  through  two  separate  methods.  Density  values  for 
the  case  above  were  calculated  at  every  grid  point  using  the 
equation  of  state.  The  alternative  method  calculated  desity 
values  at  each  station  and  allowed  the  Barnes  scheme  to 
create  a  gridded  density  field.  The  results  of  the 
alternative  method  were  very  similar  to  the  original  method. 


Table  4.1.  Mean  values  of  u  and  v  components  of true  wind 
field  and  standard  errors  for  u  and  v  components  of  sea 
level  geostrophic  winds(G)  and  terrain  following  ostrophic 
winds(  )  over  sloped  terrain  only.  Units  are  m  s 
Times  are  GMT. 


TIME/DATE 

U 

V 

U  ( G ) 

V(G) 

U(  ) 

V(  ) 

0600  21 

Dec 

2.6 

2.7 

12.6 

6.4 

10.5 

6.7 

0900  21 

Dec 

2.3 

2.8 

13.3 

7.0 

11.1 

7.4 

1200  21 

Dec 

2.  8 

2.5 

15.0 

7.7 

9.  5 

9.0 

1500  21 

Dec 

2.  8 

3.0 

14.8 

7.2 

9.7 

7.8 

1800  21 

Dec 

2.5 

2.9 

14.4 

7.6 

9.  5 

8.8 

Table  4.2  Mean  values  of  u  and  v  components  of  true  wind 
field  and  standard  errors  for  u  and  v  wind  components  of  sea 
level  geostrophic  winds(G)  and  terrain  fol^owinq  qeostrophic 
winds(  )  over  all  terrain.  Units  are.  m  s-  .  Times  are 
GMT. 


TIME/DATE 

U 

V 

U  ( G ) 

V(G) 

U(  ) 

V(  ) 

0600  21 

Dec 

3.7 

2.4 

9.5 

9.  5 

8.2 

9.2 

0900  21 

Dec 

3.6 

2.  5 

10.0 

10.4 

8.  6 

10.1 

1200  21 

Dec 

4.0 

2.6 

11.0 

11.3 

8.5 

11.6 

1500  21 

Dec 

3.6 

3.  2 

11.0 

9.8 

8.4 

9.4 

1800  21 

Dec 

3.  0 

3.  5 

10.8 

10.  6 

8.0 

11.6 

5.  Summary  and  Conclusions 

Analysis  of  surface  data  over  the  mid-Atlantic  states 
during  December  21-22,  1983,  revealed  several 
characteristics  associated  with  the  development  of  a  coastal 
front.  A  shallow  dome  of  cold  air  wedged  against  the 
southern  Appalachian  -  resulted  in  a  persistent  anticyclonic 
bulge  east  of  the  mountain  chain.  Light 

northerly ( ageostrophic )  flow  within  the  bulge  resulted  in 
light  cold  air  advection  while  southeasterly  winds  over  the 
Atlantic  produced  strong  warm  air  advection  into  the 
coastline  region.  The  net  effect  of  these  opposing  flows 
was  the  enhancement  of  the  temperature  gradient  along  the 
coast . 

Temperature  and  wind  profiles  of  individual  stations 
indicated  that  the  coastal  front  formed  along  the  coastline 
of  South  Carolina  yet  developed  over  the  interior  coastal 
region  of  North  Carolina  and  southern  Virginia.  These 
results  alter  the  scope  of  the  coastal  front  which  previous 
studies  indicated  to  be  limited  to  the  vicinity  of  the 
coast.  However,  these  findings  confirm  previous  results 
which  rejected  differential  friction  as  a  primary  force  in 
the  f rontogenesi s  process. 

Convergence  and  positive  relative  vorticity  patterns 
correlated  well  with  the  location  of  the  coastal  front  while 
divergence  and  negative  vorticity  prevailed  in  the  vicinity 
of  the  anticyclonic  bulge.  Observed  deformation  fields 


indicated  a  favorable  field  for  frontogenesis  along  the 
coastal  region.  On  the  contrary,  geostrophic  deformation 
fields  were  apparently  associated  with  the  anticyclonic 
bulge  and  had  no  relation  to  coastal  frontogenesis. 

Eulerian  and  Lagrangian  surface  frontogenesis  fields 
were  useful  in  confirming  frontogenesis  over  the  interior  of 
North  Carolina  and  southern  Virginia.  The  confluence  term 
of  the  Lagrangian  calculations  dominated  the  shear  term, 
however,  both  the  Lagrangian  and  Eulerian  frontogenesis 
fields  appeared  to  over-predict  the  increase  of  the 
temperature  gradient.  In  part,  the  overestimation  may  be  a 
result  of  the  sub  grid  scale  processes  which  do  not  allow 
the  objective  analysis  scheme  to  accurately  assess  the  true 
temperature  gradient  associated  with  the  front. 

Vertical  cross  sections  through  the  frontal  zone 
indicated  a  very  shallow  stable  surface  layer  capped  by  a 
deeper  layer  of  warmer  air.  At  the  transition  zone  from 
cold  to  warm  air,  the  vertical  temperature  gradient  was  as 
large  as  1  K  per  20  meters.  Three  dimensional  frontogenesis 
calculations  revealed  that  horizontal  confluence  and  shear 
terms  were  f rontogenetical  within  the  frontal  zone  while  the 
tilting  term  was  highly  frontolytical  within  the  frontal 
zone.  Illustrations  of  wind  components  parallel  and 
perpendicular  to  the  front  further  define  the  front  as  a 
transition  zone  from  northerly  flow  within  the  cold  wedge  to 
a  stronger  southerly  flow  which  ascends  over  the  dome  from 


the  Atlantic 


A  method  of  computing  terrain  level  geostrophic  winds 
was  used  to  test  the  hypothesis  that  pressure  reduction  was 
responsible  for  the  apparent  ageostrophic  flow  within  the 
cold  wedge.  Results  based  upon  the  terrain  level  wind 
faired  only  slightly  better  than  the  standard  sea  level 
geostrophic  winds  at  predicting  the  cold  wedge  flow  pattern. 

This  particular  case  study  afforded  an  opportunity  to 
view  an  unusual  case  of  coastal  f rontogenesis.  The  inland 
development  of  the  front  appears  to  be  an  exception  to  the 
standard  scenario.  Further  studies  of  coastal  front  events 
during  the  upcoming  East  Coast  Cyclone  Project  (GALE)  should 
attempt  to  obtain  as  much  data  on  the  vertical  structure  of 
the  coastal  front  zone  as  is  possible.  Because  this  feature 
is  so  shallow,  the  rawinsonde  observations  of  temperature 
and  winds  rarely  contain  more  than  two  readings  within  the 
cold  air  dome.  To  obtain  greater  detail  and  better 
analyses,  the  rawinsondes  should  be  adjusted  to  report 
several  readings  within  the  first  100  mb  of  the  atmosphere. 
This  detatiled  information  may  allow  future  research  to 
better  relate  vertical  temperature  and  wind  profiles  to  the 
evolution  of  the  coastal  front. 


Appendix  A 

Derivation  of  Local  Frontogenesis 
Equation  in  Cartesian  Coordinates 
As  defined  by  Pettersen { 1940 ) ,  the  local  frontogenesis 
of  any  conservative  property  S,  may  be  defined  as 

(1)  F=  dlVSI /dt 
For  any  vector 

(2)  IVSI  =(VS«VS ) *  5 

Using  the  chain  rule  of  differentiation,  (1)  becomes 

(3)  dlVSI  /at=(l/2)(VS-VS)“'*5*a(VS*VS)/at 

=  (1/(2  IVSI  ))*d(V S*VS)/dt 
The  rules  of  partial  differentiation  state 

(4)  a(A«B)/dt=A*dB/'dt+aA/at*B 

Applying  (4)  to  complete  the  differentiation  of  the  right 
hand  side  of  (3)  results  in 

aivsi  /at=i/(2  ivsi  ) * 2vs»  aivsi  /at 

(5)  =i/(  ivsi  )*vs*  aivsi  /at 

By  definition  of  the  total  derivative  and  assuming  the  total 
derivative  of  a  conservative  property  is  equal  to  zero 

(6)  aivsi  /at=VldS!/dt»V(-V«VS)+V(dS/dt) 

Substituting  (6)  into  the  right  hand  side  of  (5)  yields 

(7)  dlVSI  /dt=(  1/  IVSI  )*VS  [V( -V« VS ) +V( dS/dt ) ] 
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